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Abstract 

We consider the generation of the hierarchical charged lepton spectrum and anarchic 
neutrino masses and mixing angles in warped extra dimensional models with Randall- 
Sundrum metric. We have classified all possible cases giving rise to realistic spectra 
for both Dirac and Majorana neutrinos. An anarchic neutrino spectrum requires a 
convenient bulk symmetry broken by boundary conditions on both U V and IR branes. 
We have in particular considered the case of Majorana neutrinos with a continuous bulk 
symmetry. To avoid unwanted massless extra gauge bosons the 4D group should be 
empty. If the 4D coset is not vanishing it can provide a Wilson Line description of the 
neutrino Majorana mass matrix. We have studied an example based on the bulk gauge 
group Z7 (3)^ <8> U(3)j\f <g>j U(l) £ i with the Wilson Line in SO(3)tf satisfying all required 
conditions. A x 2 -^ to experimental data exhibits 95% CL region in the parameter 
space with no fine-tuning. As a consequence of the symmetries of the theory there is no 
tree-level induced lepton flavor violation and so one-loop processes are consistent with 
experimental data for KK-modes about a few TeV. The model is easily generalizable 
to models with IR deformed metrics with similar conclusions. 



1 Introduction 



A five- dimensional (5D) spacetime with a single warped extra dimension (WED) PQ, and two 
branes localized respectively in the ultra-violet (UV) and infra-red (IR) regions, is an exten- 
sively explored alternative to supersymmetry as a possible solution to the Standard Model 
(SM) gauge hierarchy problem, provided that the Higgs field is sufficiently localized towards 
the IR boundary. Moreover this theory has, by means of the AdS/CFT correspondence, a 
description in terms of a dual four- dimensional (4D) strongly coupled theory by which 5D 
fields localized towards the IR boundary correspond to composite states in the dual theory 
while fields localized towards the UV brane correspond to elementary ones. While both, 
supersymmetric and WED theories, can accommodate a Higgs boson with a mass around 
125 GeV, as that found by the recent LHC Higgs searches, its couplings to the different SM 
fields depend to a large extent on the different models. The comparison with experimental 
results must however wait until more accurate data on the different Higgs decay channels 
become available. 

On the other hand WED theories can successfully accommodate a solution to the flavor 
problem if all SM fields propagate in the bulk of the fifth dimension with light (heavy) fermion 
profiles leaning towards the UV (IR) boundary which are then mostly elementary (composite) 
states. This program has been very successfully applied to the quark sector in WED theories 
both for an AdS metric [21 [3], as in the original Randall-Sundrum (RS) model, and for 
asymptotically AdS (or IR deformed) metrics [IH7]. The reason of this success is that both 
the quark spectrum (in the up and down sectors) and the CKM mixing angles are hierarchical, 
a situation which can be readily described for order unity 5D Yukawa couplings if the various 
quark flavors are differently localized along the fifth dimension. In fact quark localization, 
and thus its 4D Yukawa coupling (determined by its overlapping with the Higgs profile), is 
controlled by their 5D Dirac mass c such that for different values of c (mainly in the left- 
handed doublet and right-handed singlet up sectors) we obtain hierarchically different values 
of the corresponding quark masses and small mixing angles. Moreover due to their different 
localizations the quark flavors couple differently to KK modes of gauge bosons (e.g. gluinos) 
thus generating dimension-six flavor changing neutral current (FCNC) operators. While 
these flavor violations are suppressed by the so-called RS-GIM mechanism [8] the KK masses 
have to be heavier than ~ 20 TeV to suppress FCNC processes in the RS theory [Bill] while 
this bound can be lowered to 4-5 TeV in the case of IR deformed metrics [SJE]- 

While most of the flavor literature in WED has been devoted to the quark sector, the 
lepton sector has also been extensively explored [T0l - f29] . In fact the lepton sector is quali- 
tatively different from the quark sector as the charged lepton spectrum is hierarchical while 
the neutrino spectrum and PMNS mixing angles are mostly anarchical [30]. So while differ- 
ent lepton localization is still useful to describe the charged lepton spectrum it is not that 
much in order to describe the spectrum and mixing angles in the neutrino sector. Another 
characteristic feature of the lepton sector is the neutrino nature, i.e. Dirac versus Majorana, 
which is obviously related to the possibility of lepton number violation in the bulk and/or 
the branes of the 5D theory. In fact the presence (or absence) of lepton number violation is 
mostly related to the existence of (possibly gauged because of the AdS/CFT correspondence) 
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symmetries in the bulk and the boundaries and will be the main feature of WED theories 
aiming to describe the theory of leptons. Moreover the different lepton flavor localization 
generates, by the tree-level exchange of KK modes of electroweak gauge bosons, lepton flavor 
violation (LFV) as dimension-six operators corresponding to processes as /i — > 3e and \i — e 
conversion, or in loop-diagrams from exchanged KK modes of charged leptons, neutrinos 
and the Higgs boson, in processes as fi — > e-f. Actually there is a tension between LFV 
in tree-level processes, which puts lower bounds on the Yukawa elements (the larger the 
Yukawas the more localized towards the UV brane the leptons and the more effective the 
RS-GIM mechanism), and the one-loop processes which put upper bounds on the Yukawa 
entries (the smaller the Yukawas the smaller the one-loop results as they are proportional to 
chirality- flipping mass insertions) |18] . 

In this paper we will construct a WED theory of leptons where the spectrum of charged 
leptons is hierarchical while the spectrum and mixing angles in the neutrino sector is anarchi- 
cal. We will separately review the cases where 4D neutrinos are Dirac and Majorana fermions 
and we will consider, and independently classify, all possibilities with special emphasis in 
those cases where fermion localization can successfully lead to realistic neutrino spectra. 
The former case (Dirac neutrinos) is similar to the description of the quark sector while 
the case of Majorana neutrinos is done by computing the coefficient of the dimension-five 
Weinberg operator obtained by integrating out the right-handed neutrinos in the 5D diagram 
H£i — > H£j where right-handed neutrinos are exchanged in the s-channel. This integration 
is the 5D equivalent of the seesaw mechanism in 5D theories. The general result (both for 
Dirac and Majorana neutrinos) is that if the WED theory solves the Mp/TeV hierarchy, a 
realistic anarchic neutrino sector strongly motivates the presence of a 5D symmetry which 
ensures that c\ = q and c^- = cj^ for Vz,j, where c\ and cj^- are the 5D Dirac masses for 
lepton doublets £i and right-handed neutrinos A/}. Moreover we have found for the Majorana 
case that if lepton number is broken in the bulk of the fifth dimension (whether or not it is 
broken in the UV and/or IR brane) it is not possible to accommodate the spectrum of light 
neutrinos and charged leptons simultaneously. Otherwise if lepton number is conserved in 
the bulk a realistic neutrino sector can be obtained if the theory is such that lepton number 
is violated in the UV brane, as originally proposed in [T4"] . 

The remainder of this paper is devoted to proposing a symmetry which ensures the 
anarchic structure for the neutrino sector. We will work out the case of a continuous gauge 
symmetry G in the bulk broken to the subgroup H Q (Hi) on the UV (IR) brane by boundary 
conditions. This symmetry guarantees the required structure on the 5D Dirac masses as 
well as the required Yukawa matrices for the different sectors, while in the 4D theory no 
massless gauge zero mode remains. Furthermore for the case of Majorana neutrinos we have 
proposed a Wilson Line (WL) model as the origin of the Upmns matrix which is generated 
by a background value along the coset K = K D Ki where Ki is the coset G/Hi. Using 
anarchic values of the 5D Yukawas we have made a fit to the parameters which control the 
WL and shown the 95% CL regions in the WL parameters which exhibit no fine-tuning 
in the determination of the anarchic mixing angles and neutrino masses. Our model can 
incorporate both a regular hierarchy and a inverted hierarchy of neutrino masses. Finally we 
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have considered in our model the issue of LFV. As a matter of fact, due to the structure of 
the Yukawas and 5D masses imposed by the symmetry, there is no contribution to tree-level 
processes as \x — > 3e and therefore there is no tension between tree-level and loop induced 
processes as \i — > ej. In fact the latter can be comfortably below experimental bounds for 
KK masses around 2-3 TeV and the IR Yukawa couplings of 0(1). 

The plan of this article is as follows. In Sec. |2]we have considered the case of conserved 
lepton number, i.e. Dirac neutrinos. In Sec. [3] the case of Majorana neutrinos is worked out 
in detail. The 5D propagator for right-handed neutrinos at zero-momentum is computed and 
a general expression for the coefficient of the 4D Weinberg operator explicitly given in the 
presence of lepton number violation in the bulk and in the branes for an arbitrary number of 
right handed neutrinos. In Sec. 13. II this general result where lepton number is violated in the 
bulk (and the UV and IR branes) of the fifth dimension is applied for simplicity to the case 
of only one generation of 5D right-handed neutrinos as (unlike in the 4D theories) the three 
left-handed neutrinos can receive Majorana masses by the 5D seesaw mechanism. However, 
as we have previously noticed, this case is unrealistic if one wants to describe the neutrino 
spectrum without introducing new scales, as the neutrino masses turn out to be too small. 
In Sec. 13.21 we have considered the case of an arbitrary number of right-handed neutrinos 
where lepton number is conserved in the bulk, and violated on the IR and UV branes. In this 
case the leading contribution from the previous case vanishes and the subleading correction 
can provide realistic values of the neutrino spectrum without any fine-tuning. As noticed 
above, describing the anarchic structure of the neutrino data requires a symmetry. In Sec. H] 
we have worked out a particular gauge symmetry in the bulk, consistent with an anarchic 
structure for the neutrino spectrum and mixing angles. The bulk gauge symmetry is broken 
on the branes by boundary conditions with no zero mode in the 4D theory and such that 
the matrix diagonalizing the neutrino mass matrix, Upmns, is determined by a Wilson Line 
along the coset of the broken symmetry. In Sec. [5] we have examined LFV in the present 
theory and shown that no tree-level LFV processes are generated as a consequence of the 
underlying symmetry while loop level induced processes, such as fi — > ej, are shown to give 
mild constraints on KK masses and the anarchic 5D Yukawa couplings. Finally Sec. M is 
devoted to our conclusions. We present in App. [A] details on the calculation of the general 
5D propagator of a right-handed neutrino in the presence of lepton number violation in the 
bulk and the IR and UV branes. In App. [B] we provide details of the 5D right-handed 
neutrino propagator for an arbitrary number of right-handed neutrinos, for lepton number 
conserved in the bulk and violated in the UV and IR branes. 

2 Dirac neutrinos 

Let us denote the background metric by 



ds 2 = e 2A r] flu dx fJ, dx u + dy 2 , rj^ = diag(— 1, +1, +1, +1) . 



(2.1) 



We will consider leptons as 5D Dirac fermions propagating in the bulk, with 

£i(x,y), £i(x,y), Mi(x,y) 



(2.2) 
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denoting lepton doublets, singlets and right-handed (RH) neutrinos respectively. Here we 
use % to label the three generations of ii and Si, and the (a priori arbitrary) number of copies 
of Mi- The kinetic Lagrangian for an arbitrary metric [S] reads 



Adn= [^ A {^lHl + i4>rHr) +e~ 4A (-^L + 2A> L ^ + h.c.)] 

(2.3) 



4>=e,Af,£ 

and parametrize the 5D mass Lagrangian as 

£ma SS = - Yl [ d V e ~ 4A ^ M (v) Ml + h.c. (2.4) 

tp=e,Af,s 

with constants = (— q, c/v, c^) and M(y) a function with the dimension of mass. Although 
everything can be easily worked out for arbitrary metric A(y) and mass profile M(y), we will 
for the sake of simplicity specialize to the AdS case A(y) = ky and constant mass profile 
M(y) = k. With appropriate boundary conditions there are then zero modes with profiles 

p^-c^ky p (l-2c^)kyi _ 1 

= • N *= (1-2^ < 2 - 5 > 

for ip x = £ L , £ R , Mr while the wave functions for the opposite chiralities (£r, £l, Ml) 
vanish identically. 

We will now consider a 5D Yukawa interaction between a bulk Higgs field H(x,y) and 
the leptons in Eq. (Q as| 

Ant = fdy e- 4ky £(x, y) Y N (y) M(x, y) ■ H(x, y) + h.c. (2.6) 



where H = and Y/y is a matrix of 5D Yukawa couplings that can have both bulk and 

brane contributions 

Y M (y) =Y$ + Y°5(y) + Y^5(y - Vl ) . (2.7) 
The Higgs boson zero mode profile is given by 

p aky p 2(a— 

= N,^^—^. (2.8) 

where a > 2 in the RS metric to solve the hierarchy problem. After integration over y one 
can write the neutrino mass matrix as 



my, 



v f yi s t„_„i_J 



(yjv)y e («-4-^)^ (2.9) 



1 We are using a notation where the left-handed components of the 4D SU(2)l doublets are described by 
£l while the right-handed components of the 4D neutrino and lepton singlets are described by Nr and £r, 
respectively. 

2 Of course there is a similar Yukawa interaction with charged leptons obtained from Eq. (|2.6[) by Yj^(y) — > 
Ye (y), J\f(x, y) — > £(x,y) and H(x,y) — > H(x,y), by which zero modes of charged leptons get a mass after 
electroweak symmetry breaking. 
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We will be using the fact that LH leptons should typically be leaning towards the UV 
brane in order to satisfy electroweak precision constraints, which implies c\> \. We will 
also make the reasonable assumption that 

CL ^> ~\~ Cj^f 

as a is constrained to be a > 2 as we stated above. We will compute the neutrino mass 
matrix (12.91) separately for the cases of bulk and brane Yukawa couplings. 
The neutrino mass matrix for bulk Yukawa coupling is given by 

where 

e = e~ kyi ~ 10 -15 (2.11) 
in order to solve the hierarchy problem between k ~ Mp and the TeV scale, and 



2(a- 1)11 -2c\\ |1 -2c^\ 



Yij = ~ ——j (Y£h (2-12) 

a c £ ^Af 



To avoid suppressing the r mass too much we need e^" 1 / 2 > m T /m t ~ 10~ 2 (leading to 
eg < 0.63) and therefore the case < 1/2 in Eq. (12.101) is unrealistic. On the other hand 
the case cj\f > 1 /2 in Eq. (12.101) can easily describe the spectrum of neutrino masses provided 
that Cfj > 5/6. However the neutrino spectrum predicted by Eq. (I2.10p should be hierarchical 
unless some bulk symmetry (see section H|) forces c\ = q and c]y = for Vi, which would 
render it anarchic as described by experimental data. 

In the case of an IR brane Yukawa couplings = Y^-5(y — y{) the result in Eq. (I2.10p 
applies, except that 



Y v = y/2(a - - 2^1 |1 - 2c^| (Ytf)ij , (2.13) 
while for the case of a UV Yukawa coupling Y^r = Y$5(y) we get 



e 



a-l 



(4f > 1/2) 



m» =vY 13 { , 1/9 j 7 ' (2.14) 



where 



Yy = J2(a-l)\l- 2^1 |1 - 2c^| (Y#)y , (2.15) 



At first this result seems promising as it is independent of the c\ and hence allows for large 
mixing angles El However, even for cjv > 1/2 this gives too small neutrino masses unless we 
consider e > 10~ 12 , which implies k < 10 15 GeV <C Mp, in which case we do not solve the 



3 It was previously noted that for cyv + q > a the UV Yukawa couplings are naturally dominating over 
the IR ones, leading to flavor blindness of the 4d Yukawa couplings [22] . 
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grand hierarchy problem (as in little RS models [32]). As we only focus on theories solving 
the Planck/TeV hierarchy we will disregard this class of scenarios. 

To summarize the case of Dirac neutrinos, its mass matrix and anarchic mixing angles can 
be realistically described by different localization of leptonic fields in the bulk of AdS space 
provided there is a symmetry imposing c\ = q and c\f = c^ for Vz and both ci, cm > 1/2. 
We will come back to the issue of the bulk symmetry later on in this paper. 



3 Majorana neutrinos 

The kinetic Lagrangian for one generation of RH neutrinos is given by Eq. ( 12. 3p . It is invari- 
ant under an SU(2)m global symmetry under which the neutrino transforms as a doublet 
with components M {1) = (M L ,M R ), M {2) = {M R ,-M L ) B Only M = M (1) will couple to 
Higgs and leptons, and hence SU(2)m will be broken down to lepton number L = U(1)m 
generated by a 3 . The most general bulk mass term can be written in SU(2)m covariant form 
as M^ p ■ (Tij M^ with p a real three- vector which also breaks SU (2)m- Without loss of 
generality we will choose p = (cm, 0, cm), leading to 

Anass = J e~ 4ky M(y) (-cmMrMl + h.c. + 1 -cmW l M l - MrMr + h.c.A . (3.1) 

Clearly cm conserves L, while cm breaks it. Further breaking of L can be introduced via 
the boundary conditions (BC), to be specified below. For the sake of generality, in Eq. (13. ip 
we have allowed the bulk masses to depend on y. In the following we will only consider the 
constant case M = k, while the generic case for arbitrary M(y) and A(y) is worked out in 
App.El 

The expansion in modes reads for an arbitrary number of RH neutrinos 

m*,y) ^ = W A/t )a (*,y) \ = y r ( ^ {n) W fSSiv) \ 
M a R {x,y) ) 2- ^ M R n)a (x,y) ) ^ \ W»\x) ) 1 ' 1 

where the 4D Majorana spinor (M^ n) (x),M ( - n \x)) has the 4D Majorana mass m„. At this 
point we are only considering one generation of RH neutrinos and so the label a will be 
removed. Defining new wave functions 

f%(v) = **f$(v), (3-3) 

we can rewrite the Dirac equation as 



(m n e 2k y ± c M fc)/$ = (c M k ± d y )f MRL . (3.4) 



In order to obtain the BC, notice that variation of the action Eq. (I2.3P leads to Ml(0) 
A/l(?/i) = . Including a 4D Majorana mass term for the non-vanishing field Mr, 



C 



bd 



^0 cr Cr 
2 



MrMr + h.c.l [e" 4fcyi ^MrMr + h.c. 



(3.5) 

y=yi 



4 This is sometimes referred to as a symplectic Majorana spinor. 
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(where we are considering real dimensionless numbers nf) the BC's change to 

Af L (0) + notf R (0) = , Af L ( yi ) + mtfnfa) = . (3.6) 

Clearly generic brane masses again violate L. It is important to notice that besides rij = 
also rii = oo conserve L@. 

H (o)\ / H {0) 





Figure 1: Diagram contributing to the Weinberg operator in the effective Lagrangian. 

In the general case of several generations of RH neutrinos we can always work in a basis 
where the bulk Dirac masses are diagonal. However the Majorana masses will, in general, 
be non-diagonal except if there is a bulk symmetry implying diagonal masses. Moreover in 
this work we will always impose reality of both bulk and brane mass (as well as Yukawa) 
matrices 0. We will consider the 5D Yukawa interaction between a bulk Higgs field H(x,y), 
the leptons i l (x,y) and iV/v RH neutrinos J\f a (x, y) as in Eq. (12. 6p and will integrate out 
the whole tower of RH neutrinos as in Fig. [TJ giving rise to the 4D effective operator 



i-w — c w 



if\x)-H^\x) H T ^\x)-if\x) +h.c. (3.7) 



where we denote by H <yTl \x) and if 1 (x) the normalized 4D modes of the Higgs and leptons 
doublets. To proceed we calculate the 5D propagator at zero momentum, 



a quantity which is expected to vanish whenever L is a good symmetry of the theory. Then 
the coefficient c^, with mass dimension —1, is given by 



% = / dydy' Y}?{y)Y${y>) h^\y)h^\y>) ff{y)ff{y') G% R (y,y') (3.9) 



where h^°\y) and if (y) are given in Eqs. (12. 8 p and (12. 5ft . respectively, and the neutrino 
mass matrix by 

mH = c%v 2 . (3.10) 



5 As pointed out in Rcf. [3TJ the most general BC can again be parametrized by two SUft)^ unit vectors 
Sj. The L-conserving BC's correspond to the choices s*j = (0,0, ±1). 

6 Complex entries would introduce additional CP violating phases which we are not considering in this 
work. 

7 Since for the time being we are keeping the number of RH neutrinos Nj^ arbitrary we will label them 
with indices a, b as opposed to i,j, . . . which is used for the three generations of charged leptons. 
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In the following we will investigate the implications for lepton physics in two particularly 
interesting cases. In Sec. 13. II the general case of cm ^ will be studied, while in Sec. 13.21 we 
will consider cm = 0. 



3.1 The general case cm ^ 0. 

Let us first study the case where L is violated in the bulk (and possibly on the branes). As it 
turns out one only needs to introduce one generation of 5D RH neutrinos to give Majorana 
masses to the three left-handed neutrinos. The propagator (13.81) is then given by 

c M [cq sinh(c^Q m ) - c v cosh(cQ m )] [ci smh(c v Q M ) + c u cosh(c u Q M )} fo i i v 
c v (c ci - & v ) smh(c i/ Qi) + c u (c - ci) cosh^QiJ 

where we have defined c„ = \J ' c\ + c 2 M , q = cmR; + ca/- as well as Q m = kmin(y,y'), 
Qm = kyi — kmsix(y, y') and Q\ = ky\. Details of the calculation for Grr can be found in 
App. [A] (see also Ref. [29]). 

For a bulk Yukawa coupling Yj^(y) = the coefficient of the Weinberg operator (13. 7p 
can then be written as 



-w 



+ f l 2 3 {-cv)e 2c » + f{\c u )e a ^ + f j (-c v )e ai+a ^] (3.12) 



where we have defined 



with cii 



fo( c v) = -^c v (ci + c v ){ai + c u )(aj + c v )(ai + a.j + cq - c u ) (3.13) 
fx(c u ) = Acl(cii + aj)(ai + ci)(aj + c ) (3.14) 
/^(c^) = 2c(c -c v )(a i -c I ,)(a i -c v )(a i + aj + ci + c 1/ ) (3.15) 



_ c M (a - 1) [(cq - c„)(ci + cj - (c + c 1/ )(c i - c u )e 2cv ] 1 

C u [di + CLj){ai - C v )(CLj - C v ) 



and 



f ( 24-l)6 2c H (y*) 
The leading term will be provided by the term proportional to ^(c^), resulting in 



J:j _ 2(q - l)c M e 1 (a, + a 3 - + c 1 + c^) 
(ci + c,,) (oi + a,-) 

Notice that generically (for arbitrary values of the constants a*) the rank of the matrix cy/ 
is equal to three [r(%) = 3] and, unlike in the 4D seesaw mechanism, one does not need to 
introduce several RH neutrinos to ensure that no more than one LH neutrino is massless. 
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Using the fact that LH leptons should typically be leaning towards the UV brane {c\ > 
i) and assuming Yff = 0{l)/\fk, one can estimate from (13.1 7p and (I3.18P the order of 
magnitude of the neutrino mass matrix as 

v 2 

Assuming now that c\ = q (Vi) to describe the anarchic neutrino mixing and using the 
bound e 2ci ~ l > 10 -4 imposed from the r-mass we are led to the condition m v > 10 MeV, in 
flagrant conflict with experimental data on neutrino masses 0. In the next section we will 
work out the particular case cm = where the behaviour (I3.19P appears modified by an 
extra suppression factor and one can then overcome the above problem and decouple the 
behavior of the neutrino mass matrix from that of the charged leptons. 

For an IR localized Yukawa coupling Y^(y) = Yj^5(y — yi) and using Gr_r(?/i, ?/i) = 
cm/(ci + c u ) the coefficient of the Weinberg operator (13.71) can then be written as 

where 



so that we obtain similar conclusions to the case of a bulk Yukawa coupling. Similarly, for 
a UV localized Yukawa coupling Yj^{y) = Y$-5(y) and using Gr#(0,0) = cm/(c u — c ) one 
can write 

ij _ 2(a — l)c A fe~ 2a _ 1 



— e YiYj (3.22) 
(c„ - c ) 

where 

which gives too small values for the neutrino mass m l J . Also notice that in the case where 
Yukawa couplings are exclusively localized towards one of the branes, r(cw) = 1 and the 
comments in footnote [H] do apply. 

To summarize the case where L is broken in the bulk, i.e. cm ^ 0: i) For the case of 
Yukawas in the bulk and/or localized on the IR boundary, the predicted neutrino masses 
are in the MeV range and therefore excluded by experimental data, ii) On the contrary 
for Yukawas localized on the UV boundary the predicted neutrino masses are smaller than 
experimental data by many orders of magnitude. We will see that these problems can be 
solved if the theory conserves lepton number in the bulk but breaks it on the branes. 



8 Of course this case has another flaw. If c\ = eg (Vi) it turns out that the rank of the matrix cw is 
r{cw) — 1 [ as it can easily be checked from Eq. (|3.18j) ] and the two light neutrinos are massless. This 
problem could be fixed if the symmetry enforcing equality of all c\ is approximate and c\ ~ q. However this 
case is unrealistic anyway as we have noticed because the third generation neutrino is too heavy and we will 
disregard it. 
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3.2 The case cm = 

If the 5D theory conserves lepton number L, cm = and the bulk Majorana mass term 
vanishes. In this case considering an arbitrary number of RH Majorana neutrinos one 
obtains the RR-propagator matrix as 

G RR (y, y') = e~ c ^ ky (e-^me-* - n ) _1 e~ c " kv ' (3.24) 

where = 5 ab c^j- and nf 1 are arbitrary (symmetric, dimensionless) brane Majorana mass 
matrices appearing in Eq. (I3.5|) . 

For a bulk Yukawa coupling Yj^f (y) = Y$ we write the coefficient of the Weinberg 
operator, Eq. (13. 7p . as 

Cw = y( £ c S l£ w -n o y l Y T (3.25) 

where 



Y - Y b cfr+^-e- 1 2(a-l)(24-l) 

There are different regimes, depending on the values of the numbers cjj, q and a. For 
2 < a < Cjv + q we obtain a neutrino mass scale that is suppressed at least as m u ~ 
v 2 e 2(a-i) ^ < 10~ 34 eV and hence |^ we are lead to consider the case cj\f + q < a. In this 
case assuming that the left-handed leptons lean towards the UV boundary and c/v > one 
can simplify the expression for c\y as 

c w = -Yn^Y T (3.27) 
and the neutrino Majorana mass matrix 

< ~ ^ ^K)^ (3.28) 

contains an extra suppression factor e CAr (no) _1 e CAA with respect to the cm ^ case, Eq. (I3.19p . 
which makes f!3.28j) consistent with the charged lepton spectrum [14"] . 

For Yukawa couplings localized on the IR, Yj^{y) = Y^ 5(y — yi), or UV, Y^{y) = Y$ S(y), 
boundaries the result ( I3.27P for the coefficients hold with the corresponding respective 
definitions 



-L (3.29) 
1 — e zc « 1 



^ = y » f - J *^<g-» ,3.30) 



9 In fact it has been pointed out in Ref. [22] that this case could be potentially interesting to explain the 
anarchic structure in the neutrino sector. However, as already noted there in warped space the Majorana 
case is unrealistic due to too large suppression of the neutrino mass scale. 
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We see that the case of an IR localized Yukawa coupling gives a similar expression for the 
neutrino mass as the bulk Yukawas coupling, as in the previous cases, while a UV localized 
Yukawa coupling provides unrealistic results which translate into too small values of the 
neutrino mass. 

To summarize this section we can see that the cases of a bulk or IR localized Yukawa 
coupling matrix (Y ia ) can provide a convenient description of the Majorana neutrino mass 
matrix. However the anarchic structure of the neutrino mixing angles requires that c\ = eg 
and c°j^ = cj^f for Vz, a, which requires a (gauge) bulk symmetry as we will discuss in the 
next section. In that case the general theorem on rank of matrix product implies that 
r(cw) — min{r(riQ 1 ), r(Y)} which in turn implies that to give masses to the three LH 
neutrinos we need iV/v > 3 and moreover r(riQ 1 ) > 3 and r(Y) = 3 (i.e. full rank), as in the 
cases we will consider in Sec. H] where we will fix Nj^ = 3. 

4 Flavor from Wilson Lines 

As we have seen in the previous sections the anarchic neutrino structure requires (unlike the 
hierarchical quark structure) a symmetry in the bulk. One possibility is that this symmetry 
is gauged in 5D. According to the AdS/CFT correspondence a 5D local symmetry implies 
the existence of an exact (or spontaneously broken) global symmetry of the 4D dual theory. 
Moreover, this symmetry is expected to be rather large in order to ensure degeneracy amongst 
the various Dirac bulk masses, but needs to be broken to a sufficiently small subgroup at 
the boundaries in order to allow for nontrivial Yukawa couplings. There will thus in general 
be a nontrivial coset in which some of the fifth components of the gauge fields, A5, can 
acquire VEVs. Hence, it is natural to ask if the mixing in the lepton sector could come from 
nontrivial Wilson lines or, equivalently, nonzero VEVs for A 5 . 

If we start with a group G in the bulk and break it to subgroups H and Hi at the 
boundaries, the theory has zero modes for G H fl Hi = H and A 5 e K D Ki = K with 
Ki = G/Hi. The profiles for the latter are given by 

A 5 (y) = 9e 2ky (4.1) 

where 9 is a matrix in the above coset. Each coset is spanned by the generators 

Ki = {T e GltiTHi = 0} (4.2) 

where we denote by Q (Hi) the Lie algebra of G (Hi). This is a system of dim Hi linear 
equations defining a dimG — dim Hi linear subspace. K = K Q fl Ki is the space of possible 
zero modes for A5. 

The zero mode for A§ has no potential at tree level and all configurations are degenerate, 
but can get a VEV through radiative corrections although here we will not specify the 
possible dynamics leading to different configurations. One can transform away A 5 by a 
gauge transformation. 

ryi 

A(y) = i \ A 5 (y) (4.3) 
Jy 
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The bulk action is left unchanged because of gauge invariance, provided we also transform 
all fields charged under G. This changes the UV BC of all fields transforming under the 
gauge group. In particular one has to make the replacement 

fVi 

V>(0) ->• e iA V(0) , A = A(0) = 6 e 2ky (4.4) 

Jo 

taken in the appropriate representation. The UV boundary condition for the RH neutrinos 
in Eq. fl 3 . 6 [) correspondingly changes as 

n -»■ e- iAo n e- lA ° . (4.5) 
4.1 Choice of gauge group 

We will consider the case of three generations of right handed neutrinos. The free 5D action, 
including bulk kinetic terms, is invariant under U(3)g (g) U(3)i <g> U(3)j\f. The bulk gauge 
group G should be a subgroup of the latter. In choosing G one should take into account the 
following requirements: 

• G needs to be large enough such that the breaking G — >■ H allows for nontrivial Wilson 
lines. 

• In general, the larger G, the more the theory will be protected from flavor changing 
neutral currents (FCNC). 

• G should ensure degeneracy of the ce, but allow for non-degenerate eg. 

• G should not be so large such that the breaking G — > H leaves over unwanted zero 
modes for (4D gauge symmetries). 

Since we would like to avoid hierarchical mixing, a natural choice is to take G to include 
U(3)i and to write a Wilson line to rotate RH neutrinos, we should also have G D U (3) y 
As for the charged leptons we would like to be able to write different cgi, so there should not 
be any nonabelian transformation on the charged leptons. A common charged lepton mass 
term will leave U(3)g unbroken while different masses will break it to <S>iU(l)si C U(3)s. 
The maximal group that allows to generate different charged lepton masses is thus <S>iU(l)gi. 
To summarize the simplest bulk gauge group is 

G = U(3) e ® U(3) Af ®i U(l) £i (4.6) 

where % = 1, 2, 3 runs over the number of generations. In the following we are denoting the 
11(3)^ generators by A", (a = . . . 8) where A" are Gell-Mann matrices for a = 1 ... 8 and 
A = 1. The last factor <8>iU(l)£i can be normalized such that it is spanned by {A|, A|, A°}. 
The first and second factors are consistent with constant (in generation space) q and 

10 An even more minimal choice would be [/(3)f+jv. 
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while the last one introduces differents c £ i as it is required to describe the charged lepton 
spectrum. On the other hand bulk Yukawas and the bulk Majorana mass vanish 

Y$ = Yf = c M = . (4.7) 

On the IR boundary the group G can be broken to 

H x = ®^{\) {i+£+ MY = {A| + A, 3 + A 3 ,, Al + A| + A^, A° + A° + A",} (4.8) 

i.e. "Lepton family number" which is consistent with diagonal Yukawas 

Y X>*'.r Y^ = Y £ \5 l3 (4.9) 

and forbids a Majorana mass (n\ = 0). Finally on the UV boundary the gauge group Hq is 
largely arbitrary. To allow for a Majorana mass Uq it should not contain U(l)j^, but it can 
in general contain a subgroup of SU(3)j^. Possible choices are SO(3)x, U(l)jv C SU(3)jv or 
even completely broken U(3)j^. For definiteness we will choose here the second possibility 
with the generator being X\f. The part of G acting on the doublets and charged leptons may 
be left unbroken to avoid additional Wilson line moduli in those sectors. We thus take the 
surviving group to be 

H = U{3) t ® U(l) x i u ®, t U(l) £i (4.10) 
This leads to zero UV brane Yukawa couplings 

Y$ = Y2 = Q (4.11) 

while the Majorana mass matrix no has to fulfill 

n \^ + {X^f n = . (4.12) 

in order to be Hq invariant. One can easily check that H = Hq n H% = while 

K= {Aj- 4 ' 5 ' 6 ' 7 } (4.13) 

In this case and using the previously obtained structure of Yukawas and Majorana masses 
one can write 

c w = -Y T e lA on^e iAt, Y (4.14) 
where the Yukawa coupling is given by 

Let us finally comment that a more minimal option is G = U(3)e + j^ (g), U(l) £ i with the 
boundary groups being H = U{1)i+n <S>i U{1) £ and Hi given by (14. 8p . Still H = and K 
is spanned by {A^j^-' 6,7 }. In this case Y® is nonzero and proportional to the unit matrix 
while Eq. (14.121) holds for \} +J ^ being a generator of SU(3)i + j^. As we have seen in the 
previous sections that bulk and IR Yukawas provide similar contributions to the neutrino 
mass matrix, while UV Yukawas lead to subleading contributions, the phenomenology of 
both models should be very similar. 
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4.2 A Wilson line model for Upmns 
The most general solution to Eq. (14. 12ft is given by 

n = diag ( nj, -nj, n 3 , ) (4.16) 

with rig and rig arbitrary numbers. Complex entries will result in additional phases in the 
PMNS matrix and for simplicity we will take rig 3 to be real. For the same reason we will 
only consider the subspace 

{Aj- 5 ' 7 } C K (4.17) 

which coincides with the generators of the group 50(3)^ C SU{3)m and leads to a real WL. 
For convenience we define the ratio 

2/3 = 2-f (4.18) 
such that the inverse matrix appearing in Eq. (I4.14p is proportional to 

n^<xdiag(f (4.19) 

The Yukawa matrix is diagonal because of the symmetry H\ in Eq. (14. 8p . and can be 
paremetrized as 

y ocdiag(j/i,j/ 2 ,l), (4.20) 

where we have factored out a global constant and have then normalized its last entry to 1. 
On the other hand using the fact that the coset K D 50(3)_/v we can choose the Wilson Line 
along 5*0(3) as: 

/ -b 3 i b 2 i \ 
Ao(M = hi -hi (4.21) 
V ~M hi / 

which is a linear combination of the generators \2$,i of SU(3) which span 50(3) C SU(3), 
where b k are real parameters. We can now compute the WL e lA ° as a matrix with entries 
depending on the parameters bk which in the following we will denote as 

U{b k ) = e iAo{bk) . (4.22) 

We can now use Eq. (14.141) to write the neutrino mass matrix as 

m u (b k ,y k ) oc di&g(y 1 ,y 2 , 1) • U T (b k ) ■ diag (^-, -y, lj • U(b k ) ■ diag(j/i.2/ 2 , 1) (4-23) 

where the dot indicates matrix product. The proportionality constant is given (up to 0(1) 
numbers) by e 2ci ~ 1+2cAr v 2 /e k which, using the value of q required to fix the r mass, becomes 
~ e 2cw 10 6 eV. Thus the proportionality constant is entirely controlled by and will be 
consider clS cL free parameter of the theory. 
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Once we have determined the neutrino mass matrix (14.231) we want to fit the experimental 
data consisting of the mass engenvalues nrii(bk, Uk) and the mixing angles of the matrix which 
diagonalizes the mass matrix. We parametrize the mixing angles as 

C12C13 S12C13 
U(b k ,y k )=\ -S12C23 - C12S13S23 C12C23 - S12S13S23 C13S23 (4.24) 

S12S23 ~ c 12Sl3C 2 3 -S12S13C23 - C12S23 

where we have neglected CP violation in the leptonic sector (as there are no experimental 
data on it) and we are using the notation Sij(bk,yk) = sin Qij(bk, Vk) for ij = 13, 23 and 12. 
For the mixing angles we will use the following experimental values [33] 

{sls)exp = 0.023 ±0.004 
{s\ 2 )exp = 0.312 ±0.016 

(sl 3 )exp = 0.52 ±0.06 (4.25) 

For the neutrino mass eigenvalues, for both normal < \m 2 \ < \m$\) or inverted ( 1 777.3 1 < 
\ m i\ < \ m 2\) ordering, the experimental data require 

ml — m? Am| 



r exp — ~r — t — 0.0312 ± 0.0018 (4.26) 



12 21' ex P A 2 

\ m 3 ~ m 2\ /\m A 

where Am 2 A = (2.32+°;^) x 10" 3 eV 2 and Am| = (7.50±0.20) x 10" 5 eV 2 are the atmospheric 
and solar neutrino squared mass differences respectively. The reason we only consider here 
the ratio r is that the overall normalization of the mass matrix can easily be adjusted 
to account for the correct absolute neutrino mass scale. We now want to diagonalize the 
neutrino mass (14.231) and make a fit to the mass eigenvalues and mixing angles by means of 
the x 2 function 



Notice that there is a qualitative difference between the parameters bk and as the 
former define the Wilson Line and should arise from some bulk dynamics while the latter are 
external parameters which should be 0(1) as dictated by the anarchy assumption. Notice also 
that in the case where the Yukawa matrix is proportional to the identity (i.e. yi = y 2 = ^) 
then r = 0. In that case the best (very bad) fit corresponds to U(bk,yk) = Upmns which 
happens for (61,62,63) = (0.83,0.10,0.62) and Xmin = 276. A better fit requires departure 
from 1 of y\ and/or y 2 to cope with the experimental value of the r-parameter. 

A very simple example is the case where the Yukawa and Majorana mass matrices are 
such that yk — 1. Consider for instance the case where yk = (0.90,0.95,0.90) such that 
Xmin — for 6^ = (0.84,0.11,0.62). In this case the neutrino mass spectrum at the best fit 
value is given by ~ (0.022, 0.024, 0.055) eV which has a normal hierarchy. In Fig. [2] the 
95% CL and 99% CL regions are shown in the (61,62) plane for 63 = 63 (left panel), in the 
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Figure 2: 95% CL (inside the inner ellipse) and 99% CL (inside the outer ellipse) for 
Uk = (0.90,0.95,0.90) in the plane (bi,b 2 ) with b 3 = b 3 [left panel], (b 2 ,b 3 ) with bi = b\ 
[middle panel] and (bi, b 3 ) with b 2 = b 2 [right panel]. 



(b 2 ,b 3 ) plane for bi = b® (middle panel) and in the (&i,6 3 ) plane for b 2 = b® (right panel). 
Different plots in Fig. [2] are useful to measure the available region in the space (b\,b 2 ,b 3 ) 
which is consistent with experimental data and with fixed values of Yukawas and Majorana 
mass matrix entries. 

Another different case is provided by fixed values of bk in which case we can evaluate 
the available region in the space yk which is consistent with experimental data. We will 
now provide two simple examples yielding, respectively, a normal and inverted hierarchical 
neutrino spectrum. We first consider the case bk = 0.7 (k = 1,2,3) such that Xmin — 1-3 
for = (0.30, 0.66, —0.53). In this case the neutrino spectrum at the best fit value is given 
by nii ~ (0.004, 0.010, 0.050) eV which has a normal hierarchy but more hierarchical than 
the previous example. In Fig. [3] the 95% CL and 99% CL regions are shown in the (y 2 ,y 3 ) 
plane for y\ = (left panel), in the (2/1,2/3) plane for y 2 = y 2 (middle panel) and in the 
(z/1,2/2) plane for y 3 = y° (right panel). A second example yields an inverted hierarchical 
spectrum. In this case we are considering bk = 0.4 {k = 1,2,3) such that Xmin — 2 for 
y^ = (0.60,0.63, —5.5). In this case the neutrino spectrum at the best fit value is given by 
rrii ~ (0.057,0.058,0.031) eV which exhibits an inverted hierarchy. In Fig. H]the 95% CL 
and 99% CL regions are equally shown in the (y 2 ,y 3 ) plane for y\ = y® (left panel), in the 
(2/1? 2/3) plane for y 2 = y 2 (middle panel) and in the (2/1,2/2) plane for y 3 = y® (right panel). 



5 Lepton Flavor Violation 

Warped/composite models explaining lepton masses and mixings are usually subject to se- 
vere constraints from lepton flavor violating processes. These can be mediated at tree level 
by the exchange of KK modes of electroweak gauge bosons (// — > 3e, fi — e conversion) as well 
as via loop diagrams (/x — > ej) that also involve the KK states of charged leptons, neutrinos 
and the Higgs boson. The simplest models, either with Dirac or Majorana neutrinos, are 
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Figure 3: 95% CL (inside the inner ellipse) and 99% CL (inside the outer ellipse) for bk = 0.7 
in the plane (1/2, 2/3) with y\ = y® [left panel], (yi, y-s) with y 2 = y% [middle panel] and (y\, y 2 ) 
with y 3 = [right panel]. 

required to have KK scales in the 10 TeV region for 0(1) Yukawa couplings HHEUdS] 
Moreover there is a tension between the tree- level and one- loop induced processes [IB] : while 
tree level mediated FCNC's benefit from large 5D Yukawa couplings (allowing the zero modes 
to be more UV localized/elementary, and hence to decouple from the gauge KK modes), loop 
contributions to \i — > e'y naturally grow with the Yukawa coupling. Various authors have 
thus tried to build models that reduce lepton flavor violation (LFV) via the introduction of 
either discrete [21] or continuous lepton flavor symmetries [19] 0. The model developed in 
Sec. H] has a large symmetry and can hence naturally suppress LFV processes. In fact the 
entire composite sector is completely symmetric under lepton family number since its bulk-to 
IR breaking is given by U (3)e ® U (3)./v <8>i U(l)gi — > ®iU{l)ii +£ i + ^fi and only the elementary 
sector (UV brane boundary conditions) breaks it. Since however the Higgs field is highly 
composite the elementary (UV brane localized) Yukawa interactions are completely negli- 
gible. As a consequence the charged lepton Yukawa couplings are simultaneously diagonal 
with the bulk masses for doublets (q), singlets (eg) and RH neutrinos (c/v). On the other 
hand the UV BC for the RH neutrinos breaks this symmetry and introduces the Lepton 
mixing. This can be seen in various ways 

• The UV boundary conditions are aligned with the bulk masses and IR brane Yukawas 
(no, Ytf and c\r are all diagonal), and the breaking results from a nonzero VEV of A5, 
the Hosotani breaking of the U(3)j^ symmetry. 

• The VEV for A§ is zero, Yjsf and are diagonal, but the Majorana brane mass (the UV 
boundary conditions for the RH neutrinos) becomes non diagonal, rig 1 — > U T n^U . 

n We caution the reader that the quoted numbers have to be taken with a grain of salt. Scanning over 
the anarchic 5D Yukawa couplings typically results in a broad distribution of allowed KK scales, as has been 
shown for instance in the case of the quark flavor violation [6]. Allowing for a moderate fine tuning can 
significantly reduce the bounds. 

12 We also would like to mention the possibility to reduce the bounds in the context of soft-wall models [25] . 
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Figure 4: 95% CL (inside the inner ellipse) and 99% CL (inside the outer ellipse) for bk = 0.4 
in the plane (y 2 , 2/3) with y\ = y® [left panel], (7/1, y 3 ) with y 2 = y 2 [middle panel] and (yi, y 2 ) 
with 1/3 = 7/3 [right panel]. 

• The VEV for A 5 is zero, no and c/y are both diagonal, while becomes non-diagonal, 
Y M UY M . 

These interpretations are related by a 5D gauge transformation and are completely equiva- 
lent. 

In all three cases, C£, cl and Y £ are diagonal. As a direct consequence, the mass and 
interaction eigenstates of the charged leptons are identical, all couplings to electroweak KK 
gauge bosons preserve flavor, and there are no tree-level mediated FCNC's. 




Figure 5: Diagrams contributing to u — > cy. The solid line represents a fermion which is 
either a charged lepton (a) or a neutrino (b), and the dashed line represents a boson which 
is either a neutral Higgs or Z (a) or a charged (KK) Higgs or W (b). 

As far as the one-loop contributions to \i — > are concerned, there are two types of 
diagrams, depicted in Fig. [51 depending on whether electric charge flows along a fermion 
line (a) or a boson line (b). Diagrams of type (a) never contribute in our model since all 
couplings and masses are simultaneously diagonal (no rotations on the charged fermions are 
ever necessary). To understand the contribution from diagram (b) it is most convenient 
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to work in the mass-insertion approximation and evaluate it in the basis of diagonal no 
where all flavor violation is encoded in (third point of view above). Effectively we are 
expanding the diagram in powers of the 5D Yukawa coupling. The leading contribution from 
these diagrams have one or three Yukawa induced mass insertions |27j . however it is easy 
to convince oneself that the diagrams with one insertion always involve Ye- The leading 
diagram with three Yukawa insertions is depicted in Fig. [6] and leads to the bound [27] 

2 /3TeVy ^ oQi5 

V m KK J 

where Y* is the typical size of the 5D Yukawa (in units of k) and a ~ (9(5%). Hence realistic 
KK scales of ~ 2 — 3 TeV can easily be accommodated with Y* ~ 0.4 — 0.5. Notice that, 
due to the absence of tree level FCNC's there is no lower bound on Y*. 




Y E Y N Y N 

Figure 6: The leading diagram contributing to \i — > e'y in the mass insertion approximation. 



6 Conclusion 

Warped extra dimensional models provide an elegant theory of flavor by using the fact that a 
5D Dirac mass localizes its fermion zero mode along the extra dimension in such a way that 
different localization for different fermions (i.e. different 5D masses cf) can account for the 
experimentally observed spectrum and mixing angles. Moreover FCNC higher dimensional 
operators involving light fermions (i.e. fermions localized towards the UV brane, cj > 1/2) 
and generated by exchange of gauge boson KK-modes are protected by the so-called RS-GIM 
mechanism as massive KK-modes are leaning towards the IR brane. This mechanism which 
can account for hierarchical masses and mixing angles is a very appropriate one for describing 
the mass spectrum and mixing angles (described by the CKM matrix) in the quark sector 
if the constants Cf are different in the left-handed quark doublet c qL and right-handed up 
quark singlet c UR sectors, while they are very degenerate in the right-handed down quark 
singlet Cd R sector, which suggests some protection by a flavor symmetry in the quark sector 
similar to that proposed in the present paper in the leptonic sector. 

However in the lepton sector the situation is different as the charged lepton masses are 
hierarchical while neutrino masses and mixing angles (described by the PMNS matrix) follow 
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an anarchic pattern. The localization mechanism does not work unless a fine-tuning on 5D 
masses is done, or it is implemented in a natural way by the symmetries of the theory. 
Moreover, unlike in the quark system, the nature of neutrinos, i.e. Dirac versus Majorana, is 
not yet unveiled by experiments and both situations should be considered in model building. 

In the first part of the paper we have made a systematic review of the calculation of 
the neutrino mass matrix, both in the Dirac and Majorana cases, so that we can classify 
all possible cases which could give rise to realistic spectra and mixing angles. Similarly to 
what happens in 4D if lepton number is conserved in the bulk and in both branes neutrinos 
are Dirac fermions, while otherwise they are Majorana particles. For the case of Majorana 
neutrinos we have integrated out the 5D right-handed Majorana neutrinos in the process 
£iH — > £jH, a procedure similar to the seesaw mechanism in 4D. In all cases the Yukawa 
matrix coupling the left-handed doublets li with the right-handed neutrinos Mj can be a 
5D (bulk) one and/or localized in either brane. Similarly lepton number can be violated in 
the bulk and/or in either brane. In both cases, Dirac or Majorana, the hierarchical charged 
lepton spectrum should be implemented by a pattern in the corresponding 5D masses c\. 
Moreover a realistic spectrum and mixing angle pattern for neutrinos requires: 

• For both Dirac and Majorana neutrinos: 

— A bulk symmetry implementing that c\ = q and = independently on 

— A Yukawa matrix Y^f with non- vanishing components along the bulk and/or the 
IR brane. A UV localized Yukawa matrix alone would provide too small neutrino 
masses. 

• For Majorana neutrinos: Lepton number should not be violated in the bulk. Otherwise 
the charged lepton spectrum whould lead to a too heavy neutrino spectrum. Lepton 
number violating effects are thus dominated by those from the UV brane. 

Interestingly enough the common value for the left-handed lepton doublet constant q ~ 0.63 
agrees (within la) with the common value of the right-handed down singlets c^ R — 0.65 [6] 
suggesting an extension of the flavor symmetry to the quark sector commuting with some 
unification group. 

In the second part of the paper we have constructed a simple model leading to the above 
required pattern for 5D masses, Yukawa couplings and lepton number violation. In short 
what we need is a bulk gauge group G broken by boundary conditions to the subgroup H 
(H\) on the UV (IR) brane such that: 

• The space of zero modes for A 5 , i.e. the coset space K = K fl K\ (Ki = G/Hi), is 
non-vanishing, to allow for non-trivial Wilson lines. 

• The space of zero modes for is null, i.e. H = H n Hi — to avoid unwanted 
massless non-SM gauge bosons. 

For nontrivial K one can still gauge away A 5 leading to misaligned BC for the RH neutrino at 
the two branes. The Majorana mass matrix will then depend on the WL and lead to nontriv- 
ial mixing. A priori the background (A^) is a flat direction at tree-level (a classical modulus) 
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which will however be dynamically determined at one-loop by the Coleman- Weinberg effec- 
tive potential (the Hosotani mechanism). This will then result in a dynamical determination 
of the Majorana neutrino mass matrix. Computing the one-loop radiative corrections is to 
a large extent model dependent and it is outside the scope of the present paper. Here we 
try to stress the general features of the proposed mechanism while we postpone a study of a 
model dynamically implementing the neutrino mass matrix for further studies. 

Without trying to classify all possible models we have proposed a model with G = 
U(3) e <g> U(3) M ®i U(l) £i , H x = ®iU(l) {l+s+Ny and H = U(3) e <g> U(l) x ^ ® 4 U(l) £i which 
leads to lepton number violation on the UV brane (with a particular pattern) and diagonal 
IR localized Yukawa matrices Yv an d Y £ , although we point out that other choices could do 
a similar job. In this particular example the non-trivial coset contains SO(3)m C SU(3)n 
and, for diagonal UV localized Majorana mass matrix and unit Yukawa matrix Yjy, it can 
be identified with the three angles of the PMNS rotation. Fitting also the neutrino mass 
spectrum requires a diagonal (non-unit) matrix Y/j and the Wilson line depart from the 
PMNS angles. We have quantified the available region for the Wilson line parameters and 
found that for a 95% CL region in the fit there is no fine-tuning. 

In general the larger the bulk gauge group G the more protection against FCNC's. We 
have analyzed in the present model lepton flavor violation induced by tree-level exchange of 
KK-modes (as /x — > 3e and /x — >■ e conversion) and one-loop (as /x — > ej) processes. As a 
consequence of the fact that c £ , q and Y £ are diagonal the mass and interaction eigenstates 
for charged leptons coincide, all couplings to KK gauge bosons are flavor diagonal and there 
are no tree-level lepton flavor violating processes. On the other hand lowering the rate of the 
process /x — > e^y below the experimental bound requires an upper bound on the typical value 
of the entries of the matrix Yj^. In the absence of tree-level lepton flavor violating processes 
the usual tension disappears and it is possible to find realistic values for the loop-mediated 
processes for 0(1) Yukawas and KK-masses about a few TeV. 

Finally we would like to point out that we have performed all calculations in the present 
paper for an RS (AdS) 5D metric. However a similar study would also apply for IR deformed 
metrics (as the soft-wall class of metrics [SUM]) and we expect similar conclusions to follow. 
In particular for brane localized Majorana masses and Yukawa couplings, as in the example 
we have worked out, most of the results (except for normalization factors) only depend on 
the total warp factor, which is metric independent while gauge breaking arguments should 
hold in general metrics. 
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Appendix 



A Calculating the propagator Grr for cm ^ 

In this section we give the derivation of expression (13.111) . Using the equation of motion 
(13. 4|) and the definition of G RR 



(A,; 

^— ' rn 



n 

we can write 



(M D - d y )M M \M D + d y )G RR (y,y') 

= (M D - d y )M~ M \ A f { L\y)f R n) W) +e A J2 ^rWrHV) - M M G RR (y, y') (A.2) 

n n 

Now note that the orthogonality and completeness relations are 

dye A [f L m) f L n) + fR m) f R n) } = Smn (A.3) 
e A J2[fx ] (y)&(y')] = S xx '5(y-y'). (A.4) 

n 

Assuming that there are no zero modes one can use the completeness relations to write 

[(M D - 3 y )M M \M D + d y ) + M M ] GMV, v') = Hv ~ v') (A.5) 

Now we will make the simplifying assumption Mm(?/) = CMM(y), M^iy) = cj^M(y), where 
M(y) is an arbitrary function. In other words, the 5D Dirac and Majorana masses have the 
same y dependence. Then the general solution to Eq. (1A.5|) is 

G%\ R = a^e CvQ + f3^e~ CvQ (A.6) 



where G < (G > ) refers to the regimes y < y' (y > y'), Q(y) = JqM and c v = \/ c\ + c 2 M . 
Continuity of G RR at y' gives 

(a< - a>)e c ' W) + (/?< - /3 > )e- c " Qfe ' ) = (A.7) 

The jump condition for G' gives 

(a< - « > )e c " Q(y ' ) - (/3< - /3>)e- c ^ fe,) = — (A.8) 

This can easily be solved to yield 

G< R {y,y') = b ^e c ^- Q ^+c^e c ^- Q ^ (A.9) 



G> RR (y, y') = [b+-] ^AQ(y')-Q(y)] + (c--) ^ e MQ0/)-W)l 

2 J Cy y 2 J Cy 
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The constants b and c are determined from the BC's at and y\. Note they are functions of 
y'. One can express the solutions in terms of 



min(y,y') 



Vi 



Vl 



Q, 



M. 



M , Q M = 

max(y,j/') 

Let us consider the most general BC's, which we can write as 

(c M + d y )G< R (0) = 
( Cl M + d y )G> R ( yi ) = 



Qi 



M . 



(A.10) 



(A.li; 



In the limit Cj — > oo we recover the boundary condition where fn{yi) = (here yo = 0) 
Taking q = cd however corresponds to fL{Vi) = or in other words (Mp + d y )f R = 
Applying the BC's to (IA.9I) we find the integration constants 



13 



1 c + c v 

2 AT 

1 Co — c v 

2 A r 



; Cl + c,)e c " 01 - 



- (ci + cje 1 



c,(Qi-2Q(j/')) 



(A.12) 



[ci - c^ e 



c v (-Q 1 +2Q(y')) 



N = 2 [(c ci - c£) sinh(c i ,Qi) + c l/ (c - c x ) cosh(c^Q] 

and with these we get the final expression (13. lip 

, c A/ [co sinh^Qm) - c„ cosh(c l/ Q m )] [ci suiIi^Qm) + cosh(c^Q M )] /a i o\ 
c^ (c ci - & v ) smh(c v Qi) + c v (c - ci) cosh(c„(yi) 

To give a meaning to the variables q we use the equation of motion ( 13. 4ft to eliminate /l in 
the BC's (13.61) . This results in 



(m n e A ni + (c/v + c M n;)M + c\) /r(^) = 



(A.14) 



Applying this to the definition of the propagator fl 3 . 8 [) and using the completeness relation 
(IA.4P we get 

{(cMUi + c N )M + d y ) G RR { yi ) = (A. 15) 

Comparing this to ( 1A.11I) we see that we can identify Cj as a function of the numbers Hi of 
the boundary mass terms in (13.51) . 



c M rii + cm 



(A.16) 



l.'S 



Some of these special cases have been investigated earlier [29] 
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B Calculation of Grr for n RH neutrinos and cm = 

In the case cm = and n right-handed neutrinos we get matrix equations. We will work in 
the basis of diagonal c^. The EOM (13. 4p has one term less than in the general case, due to 
the missing cm and indices are added to emphasize the matrix structure of dfj = 5 ab c b j^ and 



nf . 



m n e A f^ R a = (c a M M±d y )f^ L a (B.l) 
The BC's are given by 

fi n)a (y l ) + nff { I ? )b = (B.2) 
Multiplying the EOM (IB.ljl by & — ^and taking the sum over n we find 

e A Y,^n(v)& C {v') = (c^M±d y )J2 R/L[)X - (B.3) 

TTL n 

n n 

and using again the completeness relation one gets 

h/ R J{y - y'W c = {c a N M ± d y )G% /L Jy, y>) (B.4) 
where we have used the definition 

m n 

Solving this first order differential equation yields for Grr 

G a R b R (y,y') = e- ca ^a ab (y') (B.6) 

With Q(y) = Jq dzM(z). In contrast to the case of cm 7^ we do not have to worry about 
jump conditions as there is no delta funcion in the equation for Grr. Multiplying the BC's 
by fR(y'Y n ^ c and taking the sum over n results in a coupled equation for Glr and Grr on 
the boundary 

G LR (y u y') + n t G R R( yi , y') = (B.7) 

So in order to find a we have to find Glr first. From the EOM ( IB. 41) we find a differential 
equation with a delta function. Therefore we have to solve it for the two cases y < y' and 
y > y' and connect the resulting function with the jump condition, found by integrating the 
EOM around y' . We can write 

G LR (y,y') = GUy,y')®(y' ~y) + Gl R (y,y')e(y - y'), (B.8) 
with the stepfunction 

„, , f 1 for y > _s 

M = 1 for y < (B-9) 
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For the two cases y ^ y' we can derive 

Gf R = e c ^PHy'). (B.10) 

From the jump condition 

Gt R > (y',y')-Gt R < (y' 1 y')=6 ab (B.ll) 
the relation for the constants 0> is found: 

P< = e -cMQ(y') + /3 >_ (B.12) 

With these expressions for Glr and Grr the BC's become 

e+CArQfeO ^ e - C AAQ(j/')+/3>@( y ' _ y ) + /3 >@^ _ ^ + ¥ -w«(»i) a = o (B.13) 

which can be solved on the two branes with y^ = and result in 



a — e 



-wQ(j/i) nie - C ArQ(!/i) _ e -wQ(j'). (B.14) 



Thus we find 

G RR (y,y') = e^ Qiy) [e-^^n^-^^ - n ) e - CAfQfo,) , (B.15) 
which in the case M(y) = ky results in Eq. (I3.24p in section 13.21 
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